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Three- and four-nucleon photodisintegration processes are quite efficiently treated by 
means of effective two-body integral equations in momentum space. We recall some as- 
_-. | pects of their derivation, present previous and most recent results obtained within this 
> | framework, and discuss general features, trends and effects observed in these investiga- 
tions: At low energies final-state interaction plays an important role. Even more pro- 
lyl ■ nounced is the effect of meson exchange currents. A considerable potential dependence 
shows up in the low-energy peak region. The different peak heights are found to be 
closely correlated with the corresponding binding energies. Above the peak region only 
the difference between potentials with or without p-wave contributions remains relevant. 
,-G In the differential cross sections the electric quadrupole contributions have to be taken 

into account. The remarkable agreement between theory and experiment in p-d radiative 
capture is achieved only when incorporating this contribution, together with most of the 
above-mentioned effects. In the final part of this report we briefly review also methods 
developed, and results achieved in three- and four- nucleon electrodisintegration. We, in 
particular, compare them with a recent access to this problem, based on the construction 
of nucleon-nucleus potentials via Marchenko inversion theory. 
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1. THREE-BODY PROCESSES 

Starting form the momentum space formulation of three-nucleon collision processes, 
we sketch the transition to the corresponding photodisintegration equations. Results 
obtained within this approach are presented, general trends, sensitivities, and correlations 
are discussed. 

1.1. Deuteron photodisintegration 

The characteristic features of the momentum-space approach are most easily described 
in the two-body case. There, the scattering amplitude is given by the on-shell restriction 
(p' 2 = p 2 ) of the T-matrix 

T(p',p) = (p'\T\p) = ^(p'\V\p), (1.1) 
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where \p) represents the incoming plane-wave state, and *■ ^{p'\ the outgoing scattering 
state. The T-operator in ( |1 . 1| ) satisfies the Lippmann-Schwinger (LS) equation 

T = V + VG T. (1.2) 

Being closely related to the physical amplitude, this operator is well- defined in momentum 
space also off the energy shell (p' 2 ^ p 2 ), and the same holds true for the momentum 
representation of flTT2]), 

T(p',p) = V(p',p) + fd 3 p" V(p',p") G (p") T(p",p) . (1.3) 



There is only a singularity in the Green function Gq(j>") = (p 2 /2/i + ie — p" 2 /2n)~ l at 
the on-shell point p" 2 = p 2 , which however is easily taken care of. A useful way of doing 
this is provided by the i^-matrix approach or, more efficiently, by the VT-matrix method. 
We mention this method since it yields a particularly convenient representation of the 
two-body input to the three- and four-body treatments discussed in the following. 

Analogously to (|1.1|), the amplitude for the photodisintegration of the deuteron \ipd) is 
given by 

M(p) = ^(plHemm = (p\(l + TG )H em \^ d ) , (1.4) 

where H em denotes the electromagnetic operator. Here we have made use of the fact that 
the momentum state (p | is mapped upon the scattering state (~' (p\ via the adjoint M0ller 
operator (1 + T Go). This representation shows that the amplitude ( |1.4| ) is composed of 
a plane wave (Born) term 

B(p) = (p\He m m (1.5) 

and a further contribution (p\TGoH em \i/jd) which, via T, takes into account the final-state 
interaction(FSI). Inserting in this expression the LS equation ( |1.2| ), we obtain for M(p) 
the integral equation 

M(p) = B(p) + Jd 3 p'V(p,p') G Q {p')M{p') . (1.6) 

Equations ( |1.3| ) and ([OJ) are evidently of the same structure. Their kernels V Go, i.e., 
their most relevant ingredients, are identical. But, the nuclear Born approximation 
V(p,p') in the inhomogeneity of (|1.3| ) is replaced in (|1.6|) by the photonuclear Born term 
( |1.5| ). In other words, with this comparatively simple replacement, any working NN colli- 
sion program based on (|1.3|) can immediately be applied to deuteron photodisintegration. 
The present derivation serves primarily to exhibit structural aspects, which are typical 
also for the corresponding three- and four-nucleon generalizations. It, moreover, provides 
a quite efficient numerical tool, characterized by considerable stability, as emphasized in 
the context of ( |1.3| ). But, of course, in the two-body case there are numerous alternative 
techniques at our disposal (for a comprehensive review see [0). In the three- or four-body 
case, however, the generalization of ( |1.6|) represents one of the most efficient standard 
techniques. 



In electric dipole (El) approximation the amplitude ( |1.4|) is given by 
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where e denotes the polarization vector of the incident photon. Replacing in the com- 
mutator the free Hamiltonian H Q = P 2 /2/i by the full Hamiltonian H = H + V, this 
expression goes over into 
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with (Ef — Ej) = E 7 being the difference of the final energy Ef = p 2 /fi and the initial 
deuteron energy Ei = E d , i.e., the energy E n of the photon. 

For local potentials, which commute with X, the transition to ( |1.8|) is only another, in 
some approaches technically more convenient way of writing the amplitude (|1.7f) . For 
the typically non-local nuclear potentials, the replacement of ( |1.7| ) by (|1.8|) leads in gen- 
eral to quite different results. But, according to Siegert's theorem f2|, it is just this 
replacement, which takes into account the effect of meson exchange currents (MEC) in 
the nucleus. Siegert's argumentation involves, of course, an approximation. But, at least 
in the deuteron case it has been shown that the results obtained in this way fully agree 
with calculations based on a Hamiltonian which includes the MEC contributions explicitly 
0. In the three- and four-body cases the validity of Siegert's theorem is usually taken for 
granted, an assumption which would need a similar verification as in the two-body case. 
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Figure 1. Role of meson exchange cur- 
rents and final-state interaction. The ex- 
perimental data are from 



Figure 2. Differential cross section 

for deuteron photodisintegration with the 
Bonn B potential in El and El + El ap- 
proximation at E 1 = 60 MeV. The exper- 
imental data are from 0||. 



Fig. [I] shows the differential cross section at •& = 90° in El approximation, obtained for 
the Paris potential by solving (|1.6| ). The full solution is almost identical with the Born 



result, i.e., the FSI is practically negligible in the energy region considered. But there is 
a clear MEC effect. Fig. ^] shows the differential cross section in El approximation, and 
with inclusion of the electric quadrupole (E2) contribution. It should be mentioned that 
these results are in complete agreement with an independent solution in Ref. ||. 

One of the main questions being asked in the following is, whether similar or quite 
different FSI, MEC or El effects occur in the corresponding three- and four-nucleon 
processes. 

1.2. Three-body equations 

In the three-body case the T-matrix for rearrangement processes a + (^7) — > (3 + (0:7) 
is given by 

Tp m , an (q[], q a ) = (qp\(^f3m\Up a \tp an )\q a ) , (1-9) 

where q a denotes the momentum of particle a relative to the bound state \ip an ) of particles 
/3 and 7. Thus, due to the various initial and final configurations (a,/3 = 1, 2 or 3), we 
have, instead of a single T-operator, a set of 3 x 3 transition operators. Correspondingly, 



the single LS operator identity (|L2| ) is replaced by the coupled set of Faddeev-type Alt- 



Grassberger-Sandhas (AGS) equations [|] 

U 0a = (1 - 5 Pa )G l + £ T 7 G U ia . (1.10) 

These equations contain in their kernel the three subsystem T-operators T 7 = T a p. In 
other words, they represent in a most condensed way the relation between two-body input 
and three-body output. Of particular importance is that, like QTT2]), they are well-defined 
in momentum space, and thus amenable to the particularly accurate momentum- space 
methods. 

As three-body operators, the Up a act on the momentum state \q a ) of particle a relative 
to the (/?7) system, and on the momentum state \p a ) within this subsystem. In mo- 
mentum representation the AGS equations, hence, are two-dimensional after partial wave 
decomposition. They were used successfully in this direct form which, however, needs a 
considerable numerical effort. From the very beginning applications of ( |1.10| ), therefore, 



made use of the most typical aspect of Faddeev-type equations, the occurrence of T 7 in 
their kernel. Representing these operators in separable form, 

•*7 = 2.^/ \Xjn) t~f,nm \X7ml 1 l^-^J 

nm 

Eq. (|1.10|) reduces to a set of effective two-body equations of LS form which, in matrix 
notation, reads 

T(q', q) = V(q', q) + J d\" V{q' q") Go{q") T{q"q) . (1.12) 



Here, the matrix elements of T are off-shell extensions of ( |1.9| ), the elements of V are 
exchange potentials 

V0 m ,an(qp,qa) = (1 ~ <W) %\{Xfim\ G o\Xa.n) \(?a) (1-13) 



and the free Green function is given by 

y0;/3m,an 0/3a^a,nm ■ V / 



After this reduction, the resulting matrix LS equation ( |1.12|) is, of course, to be anti- 



symmetrized. The analogy of ( |1.3| ) and ( 1.12 ) indicates that also the photodisintegration 



equation (|1.6|) has its counterpart in the present case. A procedure, closely related to 
the steps leading to ( |1.6| ), in fact, yields for the two-fragment photodisintegration of the 
triton, 7 + t — > n + d, the effective two-body equation 

M(q)=B(q)+Jd :i q'V(q,q')go(q')M(q'). (1.15) 

Here, M.(q) and B(q) are off-shell extensions of the full photodisintegration amplitude 
M{q) = (-\q-^ d \H em \^ t ) (1.16) 

and of the corresponding plane wave (Born) approximation 

B(q) = (q\{i; d \H em \ip t ) , (1.17) 

with q the relative momentum between the outgoing nucleon and the deuteron, and \ipt) 
the triton wave function. 

The decisive advantage of this formalism is that the appropriate construction of V and 

Go, including tests of accuracy, is well known from the purely nuclear case, and has not 

to be repeated. What remains to be done, when going over to the photonuclear equation 

1.15| ), is the consistent determination of the off-shell Born term B(q). For details we 



refer to the papers quoted in the next section. 

Let us make some remarks on the derivation of ( |1.15| ). In Eq.( |1.4| ) it is the adjoint 
M0ller operator (1 + TGq) which maps (p\ onto ( ~* ) (p|. Analogously, we have in the 
present case 

M(q $ ) = (q/sKHiha + U Pa G a )H em \iP t ) , (1.18) 

with G a = (Ef + ie — H — V a )^ x the channel resolvent of the a = (/3j) subsystem. 
Proceeding in strict analogy to the derivation of (|1.6j ), i.e., inserting here, instead of the 
LS equation, the AGS equations with T 7 being chosen according to ( |1.11| ), we end up 



with an effective two-body equation of the form (|1.15|) with the matrix elements of the 
inhomogeneity being given by 

871%) = (qp\(Xpm\ G H em\i>t) ■ (1-19) 

We, moreover, mention that consistently with the approximation ( |1 . 1 1| ) the triton wave 
function may be determined by solving the homogeneous counterpart of equation ( |1.12|) . 
For details we refer to |T0|-12 . 

Three-body photonuclear integral equations were first derived within the Faddeev frame- 
work by Barbour and Phillips fl3| . A rather transparent derivation, based on the AGS 



equations, and thus more closely related to the above argumentation, was performed by 
Gibson and Lehman |14j]. It is, in fact, their approach which has been taken over in sub- 



sequent, more sophisticated applications or extensions of momentum-space techniques in 
photonuclear few-body physics. Only in parenthesis we mention that a derivation, which 



exhibits in great detail the structural equivalences pointed out above, is found in |L5[ . 

In El approximation the electromagnetic operator H em in the above three- nucleon 
relations is given by 

H'^-lf^l-P,. (1.20) 

With inclusion of meson exchange currents (MEC) we have, according to Siegert's theo- 
rem, an expression analogous to (|1.7|) , 

M(q) = -le {E f - E { ) <-><$ if, d \ £ ±±2 e • X^ t ) . (1.21) 

i=i 

Here Pj,Xj and r\ are the momentum and position operators, and the isospin Pauli 
matrices of the three nucleons involved. As in ( |1.8| ) the difference of the final energy 
Ef = q 2 /2M + Ed and the initial triton energy E{ = E t is just the energy of the incident 
photon, (Ej — Ei) = E 1 . Let us mention that for s-wave potentials only that part of the 
El electromagnetic operator which acts between the fragments survives, while, due to 
its p-wave structure, the contribution to Hex acting within the deuteron vanishes. The 
electromagnetic operator then reduces to 



H E1 = -le (Ef - Ei) ( £±i - r Ae-R 3 , (1.22) 



where R3 denotes the relative position operator between the elementary particle 3 and 
the bound state of particles 1 and 2. Although, in general, higher partial waves are quite 



relevant in realistic potentials, ( |1.22|) remains the dominant part of H em also in such cases. 
In the following considerations of the total cross section, the El contributions are 
sufficient. A considerable influence of the electric quadrupole (E2) contribution, however, 
shows up in the differential cross section. The corresponding operator is given by 

ip 3 1 _i_ T j _ _ _ _ _ 

H E2 = ~2^T,^Y JL ^' p iK-Xi + i'X j k 1 -P j ) 1 (1.23) 

with fc 7 the momentum of the incident photon. As in the above relations MEC are 
included, according to Siegert's theorem, via the replacement 

He2 = !(£/- E t ) £ l -±r L ' ■ *i ^ ■ *i ■ ( L24 ) 



1.3. Three-nucleon photodisintegration results 

First applications of the above formalism were performed in |]TB| and [|TJ|]. However, 
only the latter calculations employed the same separable Yamaguchi potential for the 
determination of V Qq and the triton wave function in B(q). Apart from a physically 
motivated different choice of the potential parameters in the bound-state case, these 
calculations, hence, represent the first consistent treatment of the problem within a simple 
separable potential model. Let us recall in this context that the Yamaguchi potential 
Vi=o(p'iP) = x(p') ^x{p)i with xip) — (p 2 + /3 2 ) -1 , is a pure s-wave potential, of course 
with different choices of the parameters A and /3 in the triplet and singlet channels, leading 
thus to an s-wave rank-one separable T-matrix. 

It is well known that, despite the simplicity of this model, the resulting n-d cross sections 
reflect essential aspects of the three-nucleon collision problem fairly well, and the same 
holds true for the corresponding photodisintegration results of |T3| and, in particular, of 



|[14|| . The main observation made already in these early solutions of the full equations 
( |1.15| ) is an about 25 % enhancement of the cross section in the low-energy peak region 
as compared to the Born approximation. This feature appears to be independent of the 
potential employed JTU] . Figure |3| shows this enhancement for the Paris potential. In view 
of the fact that the outgoing deuteron is a loosely bound, easily polarizable system, such 
a strong FSI effect is not surprising, while in the elementary-particle NN collision it was 
found to be negligible (see Sec. 1.1). Even more pronounced (about a factor of 1.7) is the 
effect of meson exchange currents. It is also much stronger than in the two-body case, an 
observation easily referred to the higher number of nucleons involved. 
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Figure 3. Role of meson exchange currents in the full solution and in Born approximation 
for the Paris potential (j < 1 + ). 



The transition operator U 0a in the nd — > nnp break-up amplitude {q',p'\Uo a \i>d)\q) is 
obtained by putting (3 = in the AGS equations ( |1.10| ) |§. This, however, means that it 



appears as a sum over the rearrangement operators Up a , and the same holds true for the 
corresponding amplitudes in the effective two-body formulation, both of the purely nuclear 
and the corresponding photonuclear processes. For details concerning the ^t — > nnp 
formalism we refer to ||16|| . Calculations along this line have been performed exclusively 



for the Yamaguchi potential in |13| and [16J . Again the more consistent calculations of |16| 
are in better agreement with experiment, although clear discrepancies are also evident. 
This calls for a repetition of these calculations with more realistic potentials. 

1.4. Potential dependence 

One of the main questions in few-body physics is the sensitivity of relevant observables 
to the underlying two-body potentials. Due to the fact that the photodisintegration 
amplitudes contain both the three-body bound and continuum states, these processes are 
expected to be a good testing ground in this respect. It, in fact, turns out that there are 
considerable differences between the photonuclear results obtained for the Yamaguchi, 
Malfliet- Tjon, Paris, Bonn A and Bonn B potentials. 

To represent these potentials or the corresponding T-matrices in separable form, various 
techniques have been developed and tested in the purely nuclear case. In the photodisinte- 
gration problem it were essentially the P^-matrix method [|17]] and the Ernst- Shakin- Thaler 
(EST) expansion 



m 



which have been employed within the above momentum-space for- 
malism. The ly-matrix method is characterized by a particularly simple separable approx- 
imation of the T-matrix, which preserves all its main analytical and unitary properties. 
The EST expansions of the Paris, Bonn A and B potentials developed by the Graz group 
(PEST, BAEST and BBEST) p|,|g, have led to the first fully reliable realistic results 

In the following applications to photoprocess an 



m 



in the three-nucleon problem 

improved parametrization by Haidenbauer is used |23 

Fig. ^ gives a good indication for the agreement between both techniques, although 
for a somewhat limited positive-parity two-body input (j < 1 + ). The results given be- 
low for the Malfliet-Tjon potential are based on the IV-matrix method. For the Paris 
and Bonn potentials, W- matrix results have also been published [[K]]. The number of 
subsystem partial waves employed in these calculations, however, turned out to be too 
small. We therefore present, instead of these previous results, most recent fully converged 
calculations ||,!1 for the PEST, BAEST, and BBEST potentials. 
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Figure 4. Comparison of the cross 

sections obtained for the Paris potential 
(j < 1 + ) using the VT-matrix and the EST 
method. 



Figure 5. Cross sections for different po- 
tentials 



According to Fig. ||, there is a remarkable potential dependence in the peak region, 
which appears to open the possibility of testing the potentials employed. Unfortunately, 
as seen in Figs. ||||, just in this sensitive region the experimental errors do not allow for 
a real judgment. This situation definitely calls for new, more accurate measurements. 
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Figure 6. Differential cross section at 
d = 90° for 7 + 3 H— > n + d obtained for 
the PEST, BAEST and BBEST (j < 1+) 
potentials. The data are from |26|- 



Figure 7. Same as Fig. |^but with (j < 2) 
contributions in the potentials. 
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Figure 8. Same as Fig. || but for the MT 
I+III and Yamaguchi potentials. 



Figure 9. Correlation of the peak heights 
and the triton binding energies for differ- 
ent potentials. 



What has to be asked for is, of course, whether this potential dependence is correlated 
with other three-body observables. Fig. |S] shows the peak heights of the cross sections 
versus the corresponding triton binding energies. In this plot the results for Yamaguchi 
potentials with different parameters and for the Malfliet-Tjon (MT I+III) potential lie 
on one straight line, the results for the realistic Paris, Bonn A and Bonn B potentials 
on another, however parallel, line. This strict separation of the two groups of potentials 



10 



is not really a statement about "realistic" or "non-realistic", but is due to the fact that 
Yamaguchi and Malfliet-Tjon are s-wave potentials, while the other ones contain p-wave 
components. In fact, switching on or off parts of these components, one can shift the 
realistic line towards the Yamaguchi-MT line p^j25[1 . In any case, the correlation exhibited 



in Fig. H is one of the typical features of few-body physics, as known from the Phillips 



23|or Tjon |HJ plots. 

At higher energies, i.e., above the peak region, the differences between the cross sections 
for the Yamaguchi and MT potentials , and between the ones for the realistic potentials 
vanish. However, there is a clear discrepancy between both groups of potentials. It should 
be emphasized that this discrepancy did not show up in our previous calculations [|l(| , due 



to the incomplete treatment of relative partial waves mentioned above. While in these 
previous calculations all potentials seemed to favor in this energy region the measurements 
by Skopik et al. |2Sj (this is still true for Yamaguchi and MT), just the realistic ones lie 
now between these data and the ones by Kosiek et al. p6| . This is another point where 
a distinct difference between potentials with and without p-wave contributions shows up. 
Also in this respect new measurements, hence, appear desirable. 

For completeness it should be mentioned that the role of p-wave contributions in the 



two-body input was emphasized already by Fonseca and Lehman [3Jj. For an investiga- 
tion of the dependence on the ci-state probability in the deuteron, with Yamaguchi-type 
interactions in these higher partial waves, see |32 . 



1.5. Radiative capture 

Detailed balance implies that the solutions of ( |1.15|) , i.e., the amplitudes for the triton 



or He photodisintegration into two fragments, determine also the inverse processes, the 
n-d and p-d radiative capture. In the above considerations of the total cross section it 
was sufficient to take into account only the El part of H em . Fig. pi] shows, in fact, that 
the El contribution, which enters the total cross section additively, is negligibly small. 

This is by no means the case for the interference term of the big El and the small E2 
amplitude occurring in the differential cross section. The E2 contribution, hence, is quite 
relevant, and evidently needed to achieve the remarkable agreement between theory and 
experiment. Let us mention that the situation is quite different in n-d capture. There, 
due to isospin factors, the interference term is strongly suppressed, and does not lead to 
any noticeable correction of the El result. 



Figure |IT] shows the same p-d process, calculated however in plane-wave (Born) approx- 
imation. The general features, in particular the E2 effect, are rather similar to the ones 
of the full calculation. But, due to the neglect of FSI, the curves lie clearly below the 
data. 



In Fig. [12] it is shown that agreement between theory and experiment is only achieved 
when incorporating the p- and /-wave contributions. Their relevance has been pointed 
out already in the photodisintegration case. To compare the corresponding observations, 
we have to bear in mind that the energy of 12.1 MeV corresponds to photodisintegration 
above the peak region. There, we also found a lowering of the cross sections for poten- 
tials with p-wave contributions which, however, led to results between the sets of data. 
The excellent agreement between theory and experiment in the present case of radiative 
capture points at the high accuracy of the experimental data, as well as the reliability of 



11 



the calculations. 

In the photodisintegration there was practically no potential dependence at somewhat 
higher energies, and thus, as expected, the difference between the Bonn A, Bonn B and 
Paris results shown in Fig. O is also rather small. 
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Figure 10. Cross section for p-d capture 
for the PEST (j < 2) potential at E = 
12.1 MeV. The data are from i3l 



Figure 11. Same as Fig. [K] but in plane 
wave (Born) approximation. 
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Figure 13. Full calculation of p-d capture 
at E = 12.1 MeV for the PEST, BAEST, 
and BBEST potentials(j < 2). 



In our review we tried to give a consistent picture of the influence of FSI, MEC, the 
potential dependence, the effect of p-waves in the interaction, or the E2 contribution in 
H em , on the cross sections of three-nucleon photodisintegration or radiative capture. It 
should be mentioned that such questions have been asked also with respect to polarization 
observables in radiative capture []3l-3~l|. In the calculations by Fonseca and Lehman, e.g., 
a strong FSI and p-wave effect was found for the tensor analyzing power T 20 and A„ 
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Very low energy radiative capture processes are of considerable astrophysical relevance. 
The n-d radiative capture, which at such energies is almost entirely a magnetic dipole 



(Ml) transition, was studied in plane wave (Born) approximation by Friar et al. [36]. In 
these investigations the authors employed their configuration-space Faddeev calculations 
of the triton wave function, with inclusion of three-body forces and pion exchange currents. 
Various trends, e.g. the correlation between cross sections and triton binding energies, 
and its potential dependence were pointed out. 

More recently a rather detailed investigation of such processes has been performed by 
Viviani et al. f37| . Their calculations employed the very accurate three- nucleon bound- 



and continuum states, obtained in the variational pair-correlated hyperspherical method 
developed, tested, and applied over years by this group. For a review we refer to the 
invited talk by Viviani at this conference ]38| . 

2. FOUR-BODY PROCESSES 

In this section we recall the generalization of the three-body momentum space approach 
to the four-body case, and review its applications to the photoprocesses 7+ 4 He — > n+ 3 He 
or p + 3 He. Similarly to the three-body case, we study the influence of FSI, MEC, the 
potential dependence, or E2 contributions on the corresponding total and differential 
cross sections. 

2.1. Formalism 

The generalization of the AGS equations ( 1.10|) is the set of 18 x 18 operator identities 



U Z = (! ~ <W) ^Go 1 ^^,,^^^ C/g 7 G T 7 G U^ (2.1) 

for the generalized four-body transition operators UZ^. As in the AGS equations, this set 
contains in its kernel the two-body operators T 7 = Tij. But now there occur also the three- 
body operators Up a given by ( pL10| ). The indices p, a, r denote the three-body subsystems 
in (3 + 1) fragmentations or the possible (2 + 2) fragmentations. As an equation for four- 
body quantities it acts on three relative (Jacobi) variables being thus three-dimensional 
after partial wave decomposition. 



Inserting the separable expansion (|1 . 1 1| ) , its dimension is reduced, as usual, by one. 



That is, we now end up with an effective two-body equation of exactly the AGS form 



LTD ) 



W p = (1 - 5 ap )g Q l + £ T T g u Tp . (2.2) 

As in the genuine three-body case it may be, and has been used as it stands, i.e., as a 
two-dimensional relation. But again it appears more convenient to further reduce the di- 
mension by expanding also the operator T T into separable terms. This two-step reduction 
scheme, concerning both the two- and three-body subsystem operators in the four-body 
equations, provides us then with an effective two-body matrix equation 

T(q>, q) = V(q', q) + J d\" V(q") G,{q") T(q", q) (2.3) 
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of the LS form (|1.3|) or (|1.12|) , of course with a much more complicated structure of the 
potentials V and Green functions Q$. Here, q and q' denote the relative momenta between 
the clusters in the (3 + 1) or (2 + 2) channels. 

And, analogously to the steps leading in the two- and three-body cases to the pho- 
todisintegration equations flTlf ) and ( |1.15| ), we end up with their four-body counterpart 
|0| 



M(q) = B(q) + Jd 3 q'V(q,q')g (q')M(q',q) 



(2.4) 



for an off-shell extension of the (properly antisymmetrized) four-body photodisintegration 
amplitude 



M(q) = 2^(q;^ in \H em \^ IV ). 



(2.5) 



Again the effective potentials and Green functions can be taken over from the purely 
nuclear case, i.e., from Eq. (2.3), and the inhomogeneity B(q) is an off-shell extension of 
the plane-wave (Born) approximation of this process, 



B{q) = 2 {q\{^ in \H em \^ IV ) . 



(2.6) 



In what follows we show results obtained for the Yamaguchi potential and the semi- 
realistic Malfliet-Tjon (MT I+III) potential. For the treatment of the two-and three-body 
subsystem operators in the kernel of the four-body equations, the ly-matrix approxima- 
tion [|I7j] and the energy dependent pole expansion (EDPE) |41j were chosen. In the 
pure nuclear case these two approximations have led to rather accurate results. The high 
accuracy of the IV-matrix representation in the photonuclear case, moreover, was demon- 
strated in Sec. 1.4 (see, in particular, Fig. 4). The accuracy of the EDPE is, of course, 
less evident and involves, in fact, one of the uncertainties of the present four-nucleon 
momentum-space approach. 

The electromagnetic El and E2 operators, with and without Siegert replacement, are 
the ones given in Eqs. ( |1.20| ), ( |1.21| ), ( |1.23| ) and ( |1.24j ). The sums runs now, of course, 
till j = 4. 

2.2. Results 

Figure 13 shows the total 4 He(7,n) 3 He cross sections for the MT potential. Due to 
the complicated cut structure of the integral equation kernel, the full calculations could 
be performed only below the three-fragment break-up threshold. Even more pronounced 
than in the three-body case (compare Fig. [3]) is the difference between the full and the 
Born, the MEC and the non-MEC results. Note that only with inclusion of these effects 
the remarkable agreement between theory and experiment is achieved [42 . 

Fig. [15] shows the corresponding results for the Yamaguchi potential. They lie below the 



MT curves, in fact too low as compared to the data. This potential dependence resembles 
completely the one exhibited in the three-body case. Again, a stronger 4 He binding (39.1 
MeV for Yamaguchi as compared to 30.1 MeV for MT I+II) leads to a lowering of the 
cross section in the peak region. 
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According to the three-body experience, FSI contributions can be neglected at higher 
energies. There, as shown in Figs. [16| and [T7|, also the potential dependence vanishes, and 
the El contributions are small. All these observations are supported by the remarkable 
agreement with newer experimental data |4"B |. The two- and three-body experiences, on 
the other hand, indicate a strong El effect in the differential cross sections 
fact, the 4 He(7,p) 3 H differential cross sections in Figs 
similarly to the ones of Figs 
remain, however, discrepancies which, according to the experiences gained in this case, 
are probably to be attributed to the neglect of higher partial waves in the two- and three- 
body subsystems. An improvement of the EDPE, the incorporation of further multipoles 
or FSI corrections may also be considered in this context. 



And, in 
T8| and ITPI show E2 corrections 
I0|, and O. In contrast to the three-body case, there 
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Figure 14. Total cross section for 

4 He(7,n) 3 He at low energies: Full solu- 
tion with exchange currents (— — ) and 

without ( ); plane wave approximation 

with exchange currents ( ) and without 

(— • — ); the data are from |43| o, [|4[ *, 
and 131 □. 



Figure 15. Same as Fig. [L4|, but with the 
Yamaguchi potential. 



2.3. Comments and alternative approaches 

The photodisintegration of 4 He into n+ 3 He or p+ 3 H has for a long time been a rather 
controversial topic. Early data appeared consistent with the picture of a giant dipole 
resonance, a feature also supported by shell model pSfl, and resonating group (RGM) 



calculations |5(| In contrast, applications of the momentum-space integral equation ap- 
proach indicated from the very beginning a fairly flat, non- resonant behavior [[40| , [5l|1 , and 
at the same time a similar trend was observed also experimentally |43[ . 



The drastic separable approximations of the three-body transition operators in the 
kernel of the four-body equations, employed in these calculations, have been replaced in 
the meantime by the more reliable EDPE, and additional new measurements are also 
available. As seen in Fig. 113, all these developments have confirmed the flat low-energy 
behavior of the total cross section predicted in [^0[ and [|3| 



As mentioned already, the RGM calculations by Wachter et al. |5(| reproduced the 
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Figure 16. Total cross section for 

4 He(7,p) 3 H at higher energies: Plane 
wave approximation including El and E2 
(— — ). Contributions of El (— — ) and 
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Figure 17. Same as Fig. [L6|, but with the 
Yamaguchi potential. 
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Figure 19. Same as Fig. [18], but with the 
Yamaguchi potential. 



apparently non-existing giant resonance. This is not surprising. Their application of the 
RGM approach is characterized by a completely unrealistic potential chosen for purely 
technical reasons. The corresponding binding energies, e.g., differ from the experimental 
ones by almost a factor of two. To cure this shortcoming, the authors replaced the 
theoretical by the experimental values in the Siegert factor. This inconsistency allowed 
them to fit first the giant resonance behavior and then, by going back to the unrealistic 
theoretical values, also the new flat results ||52| . No need to emphasize that this high 



ambiguity reduces the predictive power of such a method. 

The complexity of the full four-body integral equations suggests to look for somewhat 
less ambitious, but none the less sufficiently reliable models. The approach suggested in 



53 1 is based on writing the full outgoing scattering state as a product 



(-) 



(<f;4>. 



in 



(-) 



{q'm 



III 



(2.7) 
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Here, the two- fragment scattering state (~> (q\ describes the motion of the outgoing nucleon 
under the influence of an optical potential, adjusted to the n- 3 He or p- 3 H scattering data. 
A further simplification consists in calculating the three- and four-body bound states 
within the integrodifferential approach (IDEA) [[54]]. The results obtained in this way for 
the total cross section are in good agreement with those of the full treatment. This model, 
moreover, allows one immediately go beyond the three-body break-up threshold (compare 
the comments on Fig. [14]). The Coulomb effect, moreover, is easily incorporated. 

Note that the determination of the optical potential between the outgoing fragments 
was based in |k| on a simple analytical expression adjusted to the data. Instead, the 



more sophisticated Marchenko inversion technique could be employed, as discussed with 
respect to electrodisintegration in Sec. 3.3. 

3. THREE- AND FOUR-NUCLEON ELECTRODISINTEGRATION 

Our main topic is the photodisintegration or radiative capture of light nuclei. The 
considerable experimental and theoretical interest devoted to the corresponding electro- 
disintegration processes, however, suggests to briefly review also methods developed, and 
results achieved in this field, and to compare them with a recent alternative access to this 
problem based on Marchenko inversion. 

3.1. Three-nucleon problem 

Early calculations of electron-induced nuclear disintegration processes, with the scat- 
tered electron detected in coincidence with the ejected proton ("coincidence cross sec- 
tion"), were performed by Griffy and Oakes |55[] using the Irving three-body wave func- 
tions or a pole approximation approach. The corresponding results (dashed or dotted 



lines) are shown in Fig. |20|. Few-body techniques were employed by Lehman to obtain 
more appropriate wave functions. He calculated in this way the two- and three-fragment 
electrodisintegration of 3 H and 3 He for the Yamaguchi potential, although without FSI 
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The Gibson-Lehman approach, i.e., the Faddeev-AGS-type treatment of photodisinte- 
gration |TJ| discussed in Sec. 1.2 was extended to 3 H and 3 He electrodisintegration by 



Heimbach et a!. |p8| . Being based on the fairly simple Yamaguchi potential, this proce- 
dure takes into account consistently the final-state interaction. Van Meijgaard and Tjon 
59f employed the MT I+III potential in unitary pole expansion (UPE), and found a 



considerable difference between the Yamaguchi and MT results. 

More realistic interactions, the Urbana and Argonne potentials, were used in the varia- 
tional Monto Carlo (VMC) |JIJ calculations by Schiavilla and Pandharipande pi], while 



Laget |62[] employed a diagrammatic expansion of the amplitudes, taking into account in 



this way at least part of FSI and MEC. 

In a more rigorous approach Ishikawa et at. p3[ calculated 3 He electrodisintegration by 



a 34-channel treatment of Faddeev-type two-dimensional equations. The results for the 
Paris- and Bonn B potentials obtained in this way are shown in Fig. [H]. 

3.2. Four-nucleon problem 

In the computationally much more demanding electrodisintegration of 4 He into p + 3 H 
or n + 3 He practically no exact calculations of the final ^{q'; ipin\ scattering state, and 
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the C2 kinematics of 163 



m- 



thus of FSI, has been performed. One therefore resorts to the approximation 



^(q'-^A^l- ^(q'\(^ 



A-l 



(3. 



traditionally used in nuclear physics. Here, ^{q'\ represents the scattering state of the 
elementary particle interacting with the composite particle (ipA-i\ via an appropriately 
chosen optical potential. A decisive criterion for the quality of this approximation is that 
it does not violate the orthogonality condition between bound and scattering states. In 
other words, the condition 



(-) 



(g'l (*a-i|*a> ^ 



(3.9) 



has to be satisfied with high accuracy. We refer in this context to the investigations by 
Schiavilla and Pandharipande |JT[] and to the extension of their above-mentioned three- 
body calculations to the four-body case |64| . The three-body approach by Laget has also 
its four-body counterpart 



|65|| . The main findings of these investigations were that in 
certain kinematical regions there is an extreme sensitivity to the nuclear forces and that 
the different models employed lead to fairly distinct results. 

3.3. Marchenko inversion approach 

The considerable effort required when treating three- and four-body disintegration pro- 
cesses in their full complexity, suggests to look more closely at approximations of the type 
(|3.8| ). Assuming this approximation to be justified, the main problem consists in con- 
structing the effective interaction between the outgoing fragments in a reliable manner. 

L.L. Howell et al. |66|,67] employed for this purpose the Marchenko inversion method. 
The bound states are calculated by means of the comparatively simple integrodifferential 
equation approach (IDEA) |j||. The, in principle, exact Marchenko treatment requires 
the knowledge of the high-energy behavior of the nucleon-nucleus phase shifts. The lack 
of data at such energies allows for various extrapolations leading to different effective 
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potentials and thus, via the corresponding scattering states ^(q\, to different electro- 
disintegration results. In other words, proceeding in this way provides an appropriate 
tool for testing the high energy behavior of the phase shifts in a fairly sensitive physical 
situation. 
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Figure 22. Coincidence cross section for the two-fragment electrodisintegration of 4 He 



as a function of the missing momentum for the kinematics (A,B,C,D) of Ref. |68 



The apparent freedom in extrapolating the phase shifts practically vanishes when im- 
posing the orthogonality condition ( |3.9| ) on the resulting wave function (~>{q'\. The solid 
line in Fig. ^ obtained in this way agrees remarkably well with the pole- approximation 
results and thus with the data. In Fig. |21] the corresponding solid line lies between the 
Paris and Bonn results of Ishikawa et al. |k| . In other words, this model provides highly 
reliable results within a comparatively simple model approach, which moreover sheds some 
light on the high-energy behavior of the nucleon-nucleus phase shifts. It, hence, leads to 
quite relevant informations concerning the construction of optical potentials. 

An example for 4 He electrodisintegration at higher energies is shown in Fig. ^ The 
characteristic dip found in the plane wave (Born) approximation is partly, but not fully, 
removed when taking into account the FSI by means of Laget 's diagrammatic technique 
|65f . The Marchenko inversion approach by Howell et al. |j(],^7|] agrees fairly well with the 



corresponding results. That is, both in the three- and in the four-body case it represents 
a rather natural, comparatively simple alternative to other methods. 

With respect to the remaining discrepancies in the dip region, which are neither cured 
by the Laget nor the Marchenko inversion procedure, it appears necessary to incorporate, 
at least approximately, the interplay of the (3 + 1) and (2 + 2) amplitudes , which 
appears quite essential in the fully coupled system of equations, (compare in this context 
the discussion in Sec. 2.1). Investigations in this direction are in progress. 
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